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Abstract. We present very accurate orbital parameters 
and mass measurements (2.4% accuracy) for the well 
known very low mass triple system Gl 866. We obtain first 
orbital elements for the short-period orbit and greatly im- 
prove the long period orbit. All three stars have masses 
close to 0.1 Mq, and the system thus provides the 
strongest constraints to date on the mass-luminosity rela- 
tion close to the brown dwarf limit. 

Key words: Stars: binaries - Stars: low mass, brown 
dwarfs - Techniques: radial velocity - Techniques: adap- 
tive optics 



1. Introduction 

The nearby multiple system Gl 866 has attracted con- 
siderable attention over the last decade, as it comprises 
some of the coolest dwarfs for which dynamic masses have 
been measured. At d=3.45 pc (Van Altena et al. 1995) the 
Gl 866 system (a=22:38:34; £=-15:18:02; eq=2000.0) is a 
very close neighbour of ours, with an M5.5V joint spectral 
type (Henry et al. 1994). Leinert et al. (1986) and Mc- 
Carthy et al. (1987) independently discovered it to be a 
binary from IR 1-D speckle observations, as well as Blazit 
et al. (1987) from visible speckle observations. Leinert et 
al. (1990) first determined the elements of the orbit (P=2.2 
years) from three years of speckle observations, and Heintz 
(1993) used astrometric measurements to improve the pe- 
riod to P=2.27 years. Leinert et al. (1990) derived an ac- 
curate (~5%) total mass for the system, which however 
was too large for the measured luminosities. For the last 
few years, and at least amongst very low mass star ob- 
servers, it has been common knowledge that the brighter 
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component of the 0.35" pair is itself a short period bi- 
nary (Gl 866 AC), though the paternity of this discovery 
is unclear and, to our knowledge, any paper claim it. This 
resolves the mass discrepancy identified by Leinert et al. 
(1990). 

In this letter we present new radial velocity and angu- 
lar separation measurements, which we use to determine 
substantially improved elements for the outer orbit, to- 
gether with first orbital elements for the inner one. We 
then proceed to briefly discuss the mass-luminosity rela- 
tion for the bottom of the main sequence, in the light of 
the derived masses. 



2. Observations and data analysis 

Radial velocity measurements were obtained with the 
ELODIE spectrograph (Baranne et al. 1996) on the 1.93m 
telescope of the Observatoire de Haute Provence (France) 
between September 1995 and July 1999. This fixed con- 
figuration dual-fiber-fed echelle spectrograph covers in a 
single exposure the 390-680 nm spectral range, at an aver- 
age resolving power of 42000. The spectra were analysed 
by numerical cross-correlation with an M4V binary (i.e. 
0/1) mask, as described by Delfosse et al. (1999a). With 
this setup Gl 866 is a clear triple-lined system, with rela- 
tive depths for the three well separated correlation peaks 
of ~7, ~6 and ~1% (Figure Q). The individual radial ve- 
locities have typical accuracies of 50-80 m.s -1 for the two 
brighter components and 200-500 m.s -1 for the fainter 
one. 

Our orbital analysis largely relics on the numerous 
speckle measurements of the Gl 866 system obtained by 
Leinert et al. (1990). We have also obtained one angu- 
lar separation measurement at the 3.6-m Canada-France- 
Hawaii Telescope (CFHT) with the Adaptive Optics Bon- 
nette adaptive system PUE'O (Rigaut et al. 1998) and the 
KIR infrared camera (Doyon et al. 1998): p = .200±.001" 
and 9 = 204.1±0.9 on July 23 rd 1998, Delfosse et al. 
(1999a) provide a detailed description of the observing 
and analysis procedure, which we don't repeat here. As 
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Fig. 1. The ELODIE correlation profile on julian day 
2550389.3191 (a well separated configuration of the three 
peaks), and its adjustment by ORBIT. 

the maximum separation between components A and C is 
^0.010" (Section 3), all separation observations to date, 
whether Speckle, Adaptive Optics or with HST (Henry 
et al. 1999), only resolve the AC photocenter from the B 
components. 

We have used the ORBIT program (Forveille et al. 
1999) to determine the elements of both orbits through 
a least square adjustment to all available observations, 
including the trigonometric parallax of 0.2943" ±0.0035" 
(Van Altena et al. 1995). ORBIT directly supports triple 
systems (as well as symetrical quadruples), as long as 
three-body effects can be neglected. We have not used the 
individual velocities but instead made use of a recent im- 
provement to ORBIT (see Forveille et al. 1999 for details) 
to directly adjust the orbits to the cross-correlation pro- 
files. This direct adjustment to the profiles significantly 
improves the accuracy of the orbital parameters, since, 
for N spectroscopic measurements, it effectively removes 
~2x(N-l) free parameters per spectroscopic component. 
This gain is particularly important for Gl 866, whose three 
correlation peaks blend for many velocity configurations, 
and whose weaker C component is sometimes only de- 
tected with a moderate signal-to-noise ratio. 

In this process we produce directly a final orbit and 
not a intermediate radial velocity and we can not make 
available they, but we publish the individual radial veloc- 
ity (available in the electronic form of this paper) obtained 
by classical Gaussian fits of the correlation dips. We plan 
to send directly (by e-mail) our profile corelation at all 
people interested to work directly with them. 

3. Orbitals elements and masses of Gl 866 

Figure || and Table |^ together summarize the orbital so- 
lution which results from this simultaneous adjustment of 
both orbits to the radial velocity, angular separation, and 
parallax data. The longitude of the periastron is given 
with the spectroscopic convention, and thus refers to the 
primary The brightest and most massive component of the 



Masses 

M A (in M Q ) 0.1216 ± 0.0029 

M B (in M ) 0.1161 ± 0.0029 

M c (in M Q ) 0.0957 ± 0.0023 

Outer Orbit (AC-B) Inner Orbit (A-"CJ 

al (in a.u.) 0.775 0.012 

a2 (in a.u.) 0.414 0.016 

orbital parallaxe 290.4 ± 2.7 mas 

Table 1. Derived physical parameters. 
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system (A) and its faintest and lightest one (C) constitute 
a close binary with a ~ 3.8 day period, orbiting with the 
third star in a 823 day orbit. With P ou t/Pi„~200, the sys- 
tem is strongly hierarchical. This justifies our use of two 
non-interacting keplerian orbits in the adjustment, which 
amounts to neglecting all second order gravitational ef- 
fects. As just an example, the expected nodal precession 
period for the parameters of Gl 866, is T prec ~ 900 years 
(Harrington (1968), Mazeh and Saham (1976)), very much 
longer than the 4 years span of the radial velocity observa- 
tions (we argue below that the system is probably copla- 
nar, so that nodal precession in particular probably has a 
small or null amplitude, in addition to this long period). 
We find no secular trend in the residuals, also vindicating 
the use of two keplerian orbits. Most orbital elements are 
essentially constrained by the spectroscopic data alone, 
except for the semi-major axis, the inclination (i) and the 
orientation of the nodal line on sky (f2), which are only 
constrained by the angular separation data and the par- 
allax. 

The inner orbit has negligible eccentricity, as expected 
from tidal circularisation for its short period. The physi- 
cal separation between A and C is only 0.03 a.u. (Table 
|l|), but, thanks to the small distance to Gl 866, their ap- 
parent separation is ~0.01" . Infrared interferometers with 
baselines longer than ~40m (such as the future VLTI) 
will be able to resolve this pair, which is also an excellent 
target for observations as an astrometric binary: the large 
mass ratio and ~2 mag V band magnitude difference com- 
bine to produce a ^±0.005" astrometric amplitude in the 
passband of the FGS astrometers on HST, well within the 




Fig. 3. Phased radial velocity curves for the outer and inner orbits of G1866. In both cases the contribution of the 
other orbit has been subtracted from the displayed measurements. 
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Table 2. Orbital elements of the outer (AC-B) and inner (A-C) orbits of Gl 866. 



range of these instruments, and the outer B component 
provides an excellent local reference. 

For the time being we cannot directly determine the 
inclination of the inner orbit i s , but it is strongly con- 
strained by the requirement that Ma+Mc derived from 
the outer orbit must match the sum of the spectro- 
scopic Mxsin 3 i s obtained from the inner orbit. This 
gives sm 3 i s =0.705±0.016, and therefore i s =117°±l or 
i s ~63°±l. One of those two determinations is very close 
to the inclination of the outer orbit (i s — i\ ~ 5°). This 
probably points to a coplanar system, in keeping with a 
general tendency of close triple systems (Fekel 1981). One 
would however need to determine Q s and to resolve the i s 
ambiguity to ascertain this. 

As for Gl 570BC (Forveille et al. 1999), the combina- 
tion of high accuracy radial velocity and angular separa- 
tion data determines enough orbital parameters to derive 
the mass of each star with very good accuracy (Table |l|), 
in spite of significant uncertainties on the elements con- 
strained by the visual data alone: the system is fortunately 
close to edge-on, so that the relatively large uncertainty 
on i does not unduly propagate to sin 3 i, and the derived 
masses are fairly accurate. We obtain relative precisions of 
2.4% for both B and the AC barycenter. The short period 
spectroscopic orbit determines the mass ratio within AC 
very well, so that the individual masses of A and C also 
have this same precision. Amongst very low mass stars, 



only four systems have more accurate mass determinations 
(the three detached M-dwarfs eclipsing binaries CM Dra, 
YY Gem and GJ 2069A; and Gl 570BC). The three com- 
ponents of Gl 866 are at least twice less massive than any 
of those. Gl 866A and B are slightly more massive than 
0.1 Mq, and Gl 866C is slightly less massive. 

4. Mass-luminosity relation 

Henry et al. (1999) determine absolute magnitudes of 
M v = 15.58±0.07 for Gl 866B and (implicitly) M v = 
15.18±0.1 for the AC pair, which is too close for HSTto re- 
solve. The relative areas of the cross-correlation "dips" in 
double/triple-lined multiple systems fortunately provides 
a welcome handle on the luminosity ratios in such unre- 
solved systems. These areas are to first order proportional 
to the relative luminosities, with weaker dependences on 
effective temperature and metallicity (Mayor 1985), and 
[Fe/H] must here be identical for the three components, 
which formed from the same interstellar gas (except in 
case of binary formation by capture, wich are not prob- 
able here). For ELODIE spectra correlated with the M4 
mask, we have unfortunately not yet calibrated the spec- 
tral type (and metallicity) dependence, and as a conse- 
quence we can only derive an approximate V-band lumi- 
nosity ratio from the raw equivalent-width ratio. We ob- 
tain Lc 71^(7=0.158, which is consistent with the ^2 mag- 
nitude difference in the V band between 0.12 Mq and 
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Fig. 4. V band mass-luminosity relation. The error bars 
without points are data from Henry and McCarthy (1993), 
Henry et al. (1999) and Torres et al. (1999). The cir- 
cles represent the two well known M eclipsing binaries, 
YY Gem and CM Dra. The triangles represent our recent 
mesurements: GJ 2069A, Gl 570B and the three compo- 
nents of Gl 866. The two curves are 5 Gyr theoretical 
isochrones Baraffe et al. (1998) for two metallicities. 



0.095 M Q 5 Gyr Baraffe et al. (1998) models. We con- 
servatively assign a factor of 1.5 uncertainty to this ratio, 
which dominates the error bars on the luminosity of the 
C component. The A component on the other hand suffi- 
ciently dominates the light of the AC pair in the V band 
that even this quite large uncertainty does not propagate 
into unduly large errors on its absolute V magnitude. We 
obtain My (A) = 15.34±0.14 and My(C) = 17.34±0.45. 

Figure || represent our present best attempt at a V- 
band mass-luminosity diagram. It shows the data points 
of Henry and McCarthy (1993), Henry et al. (1999) and 
Torres et al. (1999), as well as our recent measurements of 
Gl 866, GJ 2069Bab (Delfosse et al. 1999b), and Gl 570BC 
(Forveille et al. 1999). For this last source we obtain 
absolute magnitudes of (My = 9.05±0.05 and My = 
10.89±0.05) from the V band photometry compiled by 
Leggett (1992), the HIPPARCOS parallax, and the V- 
band magnitude difference of Henry et al. (1999). The 
GJ 2069Bab eclipsing pair is clearly subluminous for its 
mass, most likely because its metallicity is significantly 
larger than solar (Delfosse et al. 1999b). The other objects 
are quite well fitted by the 5 Gyr isochrones of Baraffe et 
al. (1998), though the three Gl 866 component are slightly 
less luminous (^0.5 mag) in the V band than the solar 
metallicity model for their mass. This discrepancy most 
likely reflects a known low level problem in the models: 
the present generation of M-dwarf atmospheric models 
is specifically suspected to lack an unidentified opacity 
source in the V passband, which at the effective temper- 
ature of late M dwarfs would contribute about ^0.5 mag 
of extra absorption in this band (Allard, private commu- 



nication). Alternatively, the observations of Gl 866 can 
be reconciled with the models if its metallicity is slightly 
larger than that of the sun. With the recent and forthcom- 
ing improvements in the accuracy of M dwarfs mass mea- 
surements, metallicity determinations will increasingly be- 
come a crucial limiting factor in accurate comparisons 
with stellar models, as has long been the case for more 
massive stars (e.g. Andersen 1991). Quantitative metallic- 
ity measurements of M dwarfs are difficult in the optical 
range (e.g. Valenti, Piskunov & Johns-Krull 1998), but K 
band spectroscopy probably has the potential to reach the 
necessary accuracy (Allard, private communication). 
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Table 3. Radial velocity measurements. To be published in electronic form only, and available at the end of this 
referee version. 



